Nanoplasmonic planar traps - a tool for engineering p-wave interactions 
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Engineering strong p-wave interactions between 
fermions is one of the challenges in modern quan- 
tum physics. Such interactions are responsible 
for a plethora of fascinating quantum phenom- 
ena such as topological quantum liquids and ex- 
otic superconductors. In this letter we propose 
to combine recent developments of nanoplasmon- 
ics with the progress in realizing laser-induced 
gauge fields. Nanoplasmonics allows for strong 
confinement leading to a geometric resonance in 
the atom-atom scattering. In combination with 
the laser-coupling of the atomic states, this is 
shown to result in the desired interaction. We 
illustrate how this scheme can be used for the sta- 
bilization of strongly correlated fractional quan- 
tum Hall states in ultracold fermionic gases. 



Recently there has been growing interest in plasmonic 
nanostructures that can be used for various applications 
in quantum optics and atomic physics [H-Q ■ Particularly 
interesting is the possibility of confining atomic motion 
over regions in space of order of nanometers, comparable 
or smaller than typical values of the atom-atom scat- 
tering length. In such a regime, atomic scattering un- 
dergoes strong modifications due to confinement induced 
resonances [a> 0]- Here we propose to use this effect to 
engineer strong and robust p-wave interactions between 
fermionic atoms in planar geometries, which overcomes 
the challenges associated with creating such interactions 
in previously proposed techniques. This opens a new 
path towards the realization of exotic fractional quan- 
tum Hall states [l(| [H| and superfiuid phases [HI, EH ■ 

A strong motivation for realizing such states are their 
intriguing topological properties which find direct appli- 
cation in topological quantum computation, protected 
quantum qubits, and protected quantum memories flij ]. 
Similarly, p-wave repulsion can stabilize low filling frac- 
tional quantum Hall states [13, EB| ■> including the Moore- 
Read state [ll| . This state has been proposed in the con- 
text of a pronounced fractional quantum Hall plateau at 
filling 5/2 [l3| . but formally it also resembles the spinless 
chiral p-wavc superfiuid state. In solid-state physics, only 
in Strontium Ruthenate, chiral {p x + ip y )-w&ve Cooper 
pairs are believed to be responsible for the observed su- 
perfluidity of electrons [l6|. In the field of quantum gases, 
strong p-wave interaction can in principle be achieved 
by using Feshbach resonances. Due to the inelastic loss 
processes, however, a strong p-wave interaction is hard 



to achieve experimentally |17| . Also, in Bose- Fermi mix- 
tures, density fluctuations of bosons can induce attractive 
p-w&ve interactions or even higher partial waves between 
the fermions [I3 - l20j . However, such proposals also run 
into difficulties due to the phase separation instability of 
Bose-Fermi mixtures and stringent constraints on tem- 
perature. 

In this letter we propose to combine two important 
concepts: strongly confined two-dimensional (2D) traps 
via nanoplasmonic fields, and strong laser induced syn- 
thetic gauge fields, as illustrated in Fig.[TJ For simplicity, 
the synthetic gauge field considered is produced through 
a minimal scheme described in Rcf. |2l| . It consists of 
a laser field coupling two internal levels of the fermionic 
atoms, and an external electric or magnetic field which 
produces a linear variation of the energy of the internal 
states throughout the sample. Preparing the system in 
the lower dressed state, the sample is effectively subjected 
to a strong synthetic magnetic field [22| . Being polarized 
in one dressed state, interactions between the fermionic 
atoms are prohibited by the Pauli principle. However, 
the external degrees of freedom provide a small coupling 
to the higher dressed state. Remarkably, this will be 
shown to result in a residual p-wave contact interaction 
between the fermions. This contribution can be enhanced 
thanks to the resonant behavior of the atom-atom scat- 
tering length in strongly confined 2D settings [9(. This 
not only allows to strengthen the interaction but also to 
explore both attractive and repulsive p- wave interactions, 
i.e. going from the physics of p-wave pairing to fractional 
quantum Hall physics. 

Experimental difficulties to provide a sufficient trans- 
verse confinement, that is on the order of the atom- atom 
scattering length, arc surmountable thanks to the new 
developments in plasmonics. The interaction of cold 
atoms with nanoplasmonic systems has attracted signif- 
icant interest recently. A notable feature of surface plas- 
mon excitations, which exist along a metal-dielectric in- 
terface, is the lack of a diffraction limit. In the context of 
atom trapping, this enables the generation of fields with 
dramatically reduced effective wavelengths compared to 
free space, and a corresponding reduction of parameters 
such as trap confinement. The properties of the plasmons 
can also be greatly engineered through the underlying 
device geometry. The interaction between Bose-Einstein 
condensates and tailoredplasmonic micro-potentials has 
recently been observed [3|, and plasmon-based trapping 
techniques for ultracold atoms with applications in quan- 
tum simulation have been proposed (if, 0, Q • 
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Figure 1. The ultracold atomic sample is tightly confined above a metallic surface by a nanoplasmonic field produced by an 
external laser (in red) pointing perpendicular to the metal surface. A second laser (in blue), shining in from the side, is used 
to generate the artificial magnetic field felt by the atoms. The mechanism for fermion-fermion contact interaction involves the 
virtual excitation of one of the atoms into the excited dressed state, as illustrated in the inset. 



As a straightforward application of the scheme, we will 
consider the case of repulsive interactions and show by 
exact diagonalization how the induced p-wave interaction 
can be used to explore different quantum Hall phases, no- 
tably going from filled Landau level (LL) physics, to the 
fractional quantum Hall regime with a v = 1/3 Laughlin 
state passing through a phase with sizable overlap 
with the Pfaffian state [llj . 



Model 

We consider a trapped ultracold gas of fermionic atoms 
with two internal states |g),|e). The single particle 
Hamiltonian _ff sp = H cxt + Hal consists of an exter- 
nal part H cxt = p 2 /(2M) + V(r) with the anisotropic 
trapping potential V(r), and an atom-laser coupling Hal 
including also the internal energies. This coupling is re- 
sponsible for a synthetic gauge field which emerges due 
to the accumulation of Berry's geometrical phase when 
an atom moves within the laser field (23|. The key to 
achieving non-vanishing phases on closed contours is to 
make the internal energies, and thus Hal, spatially de- 
pendent via a Stark or Zccman shift, such that Hal and 



Hext do not commute (see Methods). 

The laser light mixes the ground and excited state, 
giving rise to position-dependent dressed states, and 
1^2), which are the eigenstates of the atom-laser interac- 
tion. As detailed in the Methods section, increasing the 
laser strength, and thus the Rabi coupling between the 
two bare states, one can energetically favor one dressed 
manifold, say l 1 !^)- By adjusting the external trapping, 
the single particle Hamiltonian projected in this lowest 
dressed manifold can be written as the usual quantum 
Hall one. 



H22 = 



(P + A) 

2M 



Mlj 2 



.(l- v 2 )(x 2 +y 2 ) (1) 



where is the effective xy trapping frequency, A = 
hr)(y,—x)/\j_, and 77 is the strength of the synthetic 
gauge field which depends on the laser wavenumber k 
and the spatial extent of the Stark or Zceman shift w 
(see Methods). This Hamiltonian has the well-known 
LL structure, and its eigenfunctions are the Fock-Darwin 
states. Restricting ourselves to the lowest LL, the corre- 
sponding wave functions read <pf B (z) cx z l exp(— |z| 2 /A5_) 
where z — x~iy describes the atom position in the (x, y) 
plane, and = \J h/(Mcu±_). The adiabatic approxi- 
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mation requires large Rabi frequencies M7o ^> En [22| . 
where the recoil energy is Er = (k 2 X 2 L /2)(hoj±). Within 
this limit, the off-diagonal Hamiltonian elements, H12 
and H21, connecting the dressed states are neglected. 
Then, transitions to the higher dressed manifold arc fully 
suppressed. A general atomic state, x( r ) = <P~i( r ) ® 
+ ^2(7*) <8> 1^2), becomes a low-lying solution for 
(pi = and (f>2 = l fif Y> - In our approach, however, some 
amount of non-adiabaticity is crucial, as it will yield a 
finite value for ipi resulting in non-zero contact interac- 
tions. 



p-wave fermion-fermion interaction 

Now we turn to the atom-atom interactions, which 
we take as contact interactions. In terms of the bare 
fermionic states, it reads 

Va = g c ^ S(z t - Zj )(\e) |g) <e| (g| + |g) |e) (g| (e|) . (2) 

Here, g c is a number quantifying the interaction strength. 
A more precise definition will be given later. Of course, 
in the dressed basis the interaction term maintains its 
form, such that interactions remain restricted to pairs 
of one atom in l^x) and the other in | \I/ 2 ) - Thus, by 
polarizing the system in the lower dressed state | \I/ 2) , no 
interactions are present in the adiabatic limit Oo — > 00. 
Still, by making the ratio of the Rabi frequency to recoil 
energy much bigger than 1, Re = KIq/Er 1, we can 
work in a quasi-polarized regime, in which the \^\) level 
serves only as a virtual manifold. 

In this limit, the unperturbed many-body Hamiltonian 
is given by 

N 

HW=Y,Hi 2 Vi. (3) 

i=l 

where the operator V% = \^2)i (^2^ projects the ith par- 
ticle onto the low-lying Hilbert space. The off-diagonal 
terms, H 12 (^2] and H 2 i ^2) and the atom- 

atom interaction of Eq. (J5|) are taken as perturbations. 
They give second-order corrections. The effective many- 
body Hamiltonian can then be written as 

H eS = H (0) + H (l) + H (2) ^ 

rri TTi 
rr(l) = V- -"21 U 12 V 

Note that the denominator in has been set to a 

constant equaling the energy difference between dressed 
states 1^2) and \*f?i). As this is taken to be large, it is 
the dominant contribution to the energy gap. 

In a previous study of a bosonic system [22|, [24j, we 
have analyzed the influence of H 1 - 1 ^, but the many-body 



contribution H ^ has been negligible due to the bosonic 
nature of the atoms. We will in the following show that 
in the fermionic case, where H 1 ^ is the only many-body 
contribution, it becomes crucial. As illustrated in Fig. [TJ 
H^ 1 describes a process where one atom is excited from 
|^2) to the virtual manifold, where it interacts with 
an atom in l^), to then get de-excited to again. 
Importantly, acting among fcrmions, the many-body in- 
teraction term gives solely non-zero p-wave contributions. 
This is seen by using Eq. (TT3"|) to cast into p-wave 
form (cf. Ref. (H), 

ffP'a^^^^^P^ (5) 

with the relative variables, pij = —ih(d Zi — d z ) and z,j = 
Zi — Zj. The important feature of Eq. Q is that the 
effective interaction is linear in the bare one Vij, which 
allows one to change the interaction from attractive to 
repulsive. This is in contrast to second order mechanims 
like the Kohn-Luttinger (25| . 

The main question which arises at this point is whether 
the residual interaction term, Eq. ([5]), is strong enough 
to significantly modify the physics of the system. This 
becomes possible by tuning the interaction strength g c . 
It is well known that this parameter crucially depends on 
the geometry of the system. In particular, for transversal 
confinements on the order of the scattering length, and 
considering the case of attractive interaction the effective 
2D coupling is known to behave as 0] 



M V2^\ z /a 3D + log(0.918fto; z /7re) ' 

where e is the energy of the motion in the x — y plane and 
a,3D the 3D scattering length. For a value of hu> z /e = 10 3 , 
it produces a resonant behavior for values of the trans- 
verse confinement X z ~ 0.4|a3£>|. This confinement- 
induced resonance behaviour is not present in usual ex- 
periments with optical traps. There, the trapping on the 
z direction has at most been of the order of hundreds of 
nanometers, far from the resonance region. The trans- 
verse confinement lengths of X z ~ 5-10 nm needed to 
facilitate significant interactions can be achieved using 
novel plasmon-bascd trapping techniques, such as those 
investigated theoretically and experimentally in Refs. [3- 
[j. For example, it is possible to tailor a two-dimensional 
array of metallic nanosystems (such as nanoshells Q), 
which creates a near-planar trapping potential arising 
from spatial interference between an incident field and 
plasmon-enhanced near-field. The effective wavelength 
characterizing this trapping potential scales like the char- 
acteristic size of an individual nanosystem, X e g ~ r, 
even for system sizes far below the free-space wave- 
length r <C Xq. This yields a corresponding reduction 
of ~ ^/r/Ao in the trap spatial confinement compared to 
free-space techniques. 

This resonant behavior can in principle be used to pro- 
duce arbitrarily large values of g c and, importantly, al- 
lows to achieve not only large values of the coupling, but 
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Figure 2. Evolution of the ground state as a function of the in- 
teraction parameter g c = g c /(k\± ijRe) 2 - Panel (a) presents 
the overlap of the ground state with the filled LL state (solid- 
black), the fermionic Pfaffian state (dotted-red), the quasi- 
particle state over the v = 1/3 Laughlin (dashed-orange) and 
the v = 1/3 Laughlin (long-dashed- blue). Panel (b) contains 
the average angular momentum of the ground state of the 
system. 

also provides a way of producing both attractive and re- 
pulsive p-wave interactions between the fcrmions. 

Example: stabilization of the v = 1/3 Laughlin state 

As an example, we discuss the case of repulsive p-wave 
interaction. In the context of quantum Hall physics we 
have to ask whether the obtained p-wave interaction is 
capable of bringing the system from the integer quantum 
Hall regime of a non-interacting system to the fractional 
quantum Hall regime. In that case, a Laughlin-like state 
should show up as the ground state of the system 0- Note 
that in the quantum Hall regime, 77 — >■ 1, the contribution 
of i?22 reduces to a constant, as all Fock-Darwin states 
become (quasi)degenerate. Thus, to bring the system 
into the fractional quantum Hall regime, the interaction 
term must be comparable to the contribution of H^- 1 ' 
term, which breaks the rotational symmetry [24j . 

To give definite numerical predictions, we perform an 
exact diagonalization (see Methods) with a few number 



1 Let us recall the Laughlin wave function llOll at filling is, 

^Laughlin = li fe " ^ ) 1/ ^- !z|2/2 • (7) 
i<j 

which has L z = u~ 1 N(N - l)/2. 



of atoms, N = 4. The parameters of the system are taken 
as k = 10/Aj., Hflo = 100E R and 77 = 0.98. Wc discuss 
the different phases appearing as we vary the interaction 
strength, g c . For weak interactions, the ground state of 
the system has a large overlap with the analytical form of 
the filled LL state, v = 1, as depicted in Fig. [2] The angu- 
lar momentum of the ground state is found to be slightly 
larger than the analytical value, L = 6. As explained in 
Refs. [22], HH, this is due to the derivation from rotational 
symmetry. When the interaction is increased, the system 
undergoes a transition into a phase, where the ground 
state has large overlaps with the fermionic Moore-Read 
state [ll|. Even stronger interactions bring the system 
into a state which resembles the quasiparticle excitation 
of the v = 1/3 Laughlin state. At another critical value 
of g c , one finally reaches the 1/ = 1/3 Laughlin state. As 
in the case of bosons with contact interactions, this state 
has zero interaction energy, and thus for any stronger 
interaction parameter, it remains the ground state. 



Summary 

We have presented a novel mechanism to realize sizable 
p-wave interactions between fermionic atoms. The key is 
the combination of a strongly confining plasmonic field, 
which allows to explore confinement-induced resonances, 
with a simple scheme to generate a strong synthetic gauge 
field. To exemplify the potential of our approach, we have 
considered the case of repulsive p-wave interaction. We 
have shown that our proposal allows to stabilize a num- 
ber of interesting quantum Hall states, like the Pfaffian, 
and the v = 1/3 Laughlin state. In our numerical cal- 
culation we have considered a small number of atoms, 
as has become experimentally feasible recently (26l - [28| . 
but we note that the scheme should also be applicable 
to large systems. A good candidate for realizing the pro- 
posal are Ytterbium atoms due to the long-lived state 
of the clock transition. Requiring an ultratight trapping 
in a 2D geometry, our proposal shall trigger the use of 
nanoplasmonic fields as a promising technique for achiev- 
ing that goal. 
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METHODS 



Single particle Hamiltonian 



Ref. [22j. Note that for w <C x,y, we recover the symmet- 
ric gauge expression A(r) = j^(y,—x). With a conve- 
nient choice of the trapping potential, H22 can be made 
symmetric. Then, the Hamiltonian element H22 reads 



H 



22 



ip + A) 



Muj 2 



(1 -r] 2 )r 2 



(11) 



The single particle Hamiltonian reads, 



H sp — -ffext 



Hal 



(8) 



2M 2 

where lu±_ is the effective xy trapping frequency, r\ = 
(kX 2 L )/(Aw), and Aj_ = y/H/(Mu±). The recoil energy 
of the atoms is defined as Er = (k 2 X 2 L /2)(hu>±). 

Retaining up to quadratic terms, the off-diagonal 



where H c ^ t = P 2 /(2M) + V(r) with the anisotropic trap- 
ping potential V(r). Hal is the atom-laser coupling, 
which includes the internal energies. To make Hal spa- 
tially dependent, such that Hal and H ext do not com- 
mute, we perform a Stark or Zeeman shift of the internal 
energies. The strength of this shift can be characterized 
by a length scale w, which is chosen such that the energies 



Hamiltonian elements, iJ 12 

h 2 



Hi-,, explicitly read 



H 



12 



with a T 



2M 
' AM 



k x 
Aw 



ac\ 



+ a/C2 + bc 3 + 



-X±dg + A 



-1 1 



of the bare internal states read E g = — MIqx/(2w) and — X±d z + X ± hz, and b = X±d. 



iky 
Aw 2 

X±d, 

+ x 



: + X 
1 1 „ 



w 



E e = hujA + MIqx/(2w). Here, uja is the energy differ- 
ence of the bare states. In this way, preparing the system 
in the ground-state of Hal, the external part will stim- 
ulate transitions into the excited manifold of Hal- The 
probability of such transitions is controlled by the Rabi 
frequency Qq of the coupling. The laser frequency is set 
to resonance with the atomic transition. Furthermore, we 
choose the laser to be a running wave in y-direction with 
wavenumber k. Then, within the rotating-wave approxi- 
mation, the atom-laser Hamiltonian Hal can be written 
in terms of bare states |e) and |g) as [23 



(12) 

St EE 



Acting on 



-1 1 . 



H 



AL 



~2~ 

4- sin 



cos0(|e) (e| 



(e^ |e) 



h.c 



;gi) 
)]. 



a Fock-Darwin state the operators a (a') decrease (in- 
crease) the I quantum number by one, while the opera- 
tors b and ¥ change the Landau level. 

As we will be interested in the fractional quantum 
Hall regime of large synthetic magnetic field, r\ ~ I, 
it is possible to safely neglect the b and ¥ contribu- 
tions. In this limit, we have c\ = c 2 ~ %wjX\, and 
C3 = — C4 ~ 2/(iyAj_). In our numerics, we will further- 
more choose w ~ 2.5Aj_ and k = 10/Aj_, implying rj ~ I, 
guaranteeing c\ -C C3. We can then write 

H12 = "^(fi + « f ) + 0[(w/X ± )-% (13) 



(9) 



where fl = f2oyl + x 2 /uu 2 , tan# = w/x, and <j> = ky. 
Note that spontaneous emission processes are not con- 
sidered in the Hamiltonian of Eq. ([9]). This is justified 
if the two atomic states are sufficiently long-lived, as is 
the case for the 1 S'o — s- 1 Pi clock transition in Ytterbium. 
In contrast to the bosonic Ytterbium isotopes, the finite 
spin of the fermionic isotopes yields a small magnetic 
moment, which allows for a strong coupling of the clock 
states at reasonable laser power. Thus, achieving large 
Rabi frequencies, as required by our proposal, poses no 
problem for 171 < 173 Yb 

Diagonalizing Eq. © yields the dressed states, 
= c- i;G (Ce^/ 2 |g) + Se-^/ 2 |e)), |* 2 ) = 
c iG (-S*c^/ 2 |g) + Ce-W 2 |e)), wherc c = cos6 )/ 2 , 

S = sin 9/2, G = The single- particle Hamiltonian 

H sp can be expressed as a 2 x 2 matrix i?y . In the dressed 
state basis, its diagonal terms can be written as (22|, 



Exact diagonalization 

To solve the effective Hamiltonian, we perform exact 
diagonalization. Therefore, we build many-body states 
using as single particle states the Fock-Darwin states, 
\l). Then the second quantized form of H^> is 



a 2 \Y,WM V < 



ij : kl 



(14) 



ij,kl 



where anihilates an atom in ipf B (z). The matrix el- 
ement reads, V ljM = {MIq)-' 1 j\H 2 iV Hi 2 \l)\k) . Tak- 
ing into account the Pauli principle, we get, 



(2) 



(<p*h*ip k hi) - (ip*h*h k <pi) 



-(h*ip*jip k hi) + (h*<p*h k (pi) 



(15) 



H, 



2M 



U + V + e^, 



with ei = 1 and e 2 = — 1- Full expressions for the vec- 
tor potential A and the scalar potential U are given in 



where hi = [Hi 2 / (hQ Q )] tpf . With the expression for 
H\2 from Eq. (fT"3"|) . hi is directly found to be hi = 
(kX± t]Re)~ 1 (VI + 1 <Pl+i + tpi-x) ■ It is worth noting 
that due to the contact nature of the interaction, Eq. (fT4")) 
commutes with L 7 . 
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